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ABSTRACT 

We investigate the constraint on the abundance of primordial black holes (PBHs) and the 
spectral index n of primeval density perturbations given by the ionizing photon background 
at the epoch of reionization. Within the standard inflationary cosmogony, we show that the 
spectral index n of the power-law power spectrum of primeval density perturbations should be 
n <1.27. Since the universe is still optical thick at the reionization redshift z ~ 6 - 8, this 
constraint is independent of the unknown parameter of reheating temperature of the inflation. 
The ionizing photon background from the PBHs can be well approximated by a power law 
spectrum J(z^) oc v^, which is greatly different from those given by models of massive stars and 
quasars. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. INTRODUCTION 

Primordial black holes (PBHs) form from the grav- 
itational collapse when an initially dense region on 
super-horizon scales enters the horizon during the 
early radiation-dominated era (Carr & Hawking 1974, 
Carr 1975). In the inflationary universe, the abun- 
dance of PBHs is extremely sensitive to the spectral 
index n, which specifies the dependence of the power 
spectrum of primordial density fluctuations on the 
comoving wavenumber /c, oc A;". Thus, a con- 

straint on the abundance of PBHs will also yield a 
constraint on the index n. 

The simplest constraint on the abundance of PBHs 
is given by the condition ri^/j < 1, i.e., the mass den- 
sity of the PBHs should be less than the critical den- 
sity of the universe. This condition results in n < 1.5 
regardless of whether a black hole can evaporate com- 
pletely, or ends with a Planck relic (Kim & Lee 1996). 
More effective constraints can be derived by com- 
paring the particle emission from PBH evaporation 
with observed background. For instance, the cosmic 
background of 7-ray spectrum on ~ 100 Mev places 
a tightest constraint on the PBH density parameter 
^bh < (MacGibbon k Carr 1991), which yields 
the index n < 1.25 (Kim, Lee & MacGibbon 1999). 
However, these constraints depend on an additional 
parameter, the reheating temperature T^/j, or equiv- 
alently, the minimum mass of PBHs M^in, which is 
not well determined. 

Recently, the epoch of the onset of the cosmic reion- 
ization is explored with absorption spectrum of high 
redshift (5.7 < z < 6.3) SDSS quasars (Fan et al. 



2001, and references therein). The existence of the 
complete Lya Gunn-Peterson trough in the absorp- 
tion spectra of high redshift quasars indicates that the 
dark age probably ended at the redshift z ~ 6 — 8. 
The ionizing photon background in the IGM at z ~ 6 
is found to be more than 20 times lower than that 
at z ~ 3. These results motivate us to explore the 
constraints on the abundance of PBHs, and then on 
the primordial density perturbations with the low en- 
ergy photon background at the epoch of reionization. 
We will show that this constraint is basically inde- 
pendent of the uncertain parameters Trh or 
PBHs have been considered possible source of 
the cosmic reionization (Gibilisco 1996, 1997), but in 
these works, the mass function of PBHs is treated as 
free parameters. It cannot be used to constrain n. 

This Letter is organized as follows. In §2, we briefly 
present the methods of calculating the mass func- 
tion of the PBHs and the ionizing photon background 
from PBH's evaporation. The numerical results of the 
constraints on n from the observed high redshift ion- 
izing background are given in §3. §4 are discussions 
and conclusions. 

2. METHODS 

2.1. Mass function of PBHs 

For Gaussian primeval density perturbations, the 
mass function ni,h{Mbh) of PBHs at redshift z can be 
calculated by the same way as calculating the mass 
function of collapsed massive halos with the Press- 
Schechter formalism (Press & Schechter 1974). Ac- 
cordingly, the number density of regions with density 
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contrast from 6 to 6 + d6 and mass M to Af + dM at 
the time ti of the onset of radiation-dominated epoch 



is 



n(M, 6)dMdS 



(1) 
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where pi is the cosmic mean mass density at ti, M 
the mass within R, and aji the variance of the prob- 
abihty distribution function of initial Gaussian mass 
field smoothed by scale R. Since the mass function is 
normalized at pi, i.e. / / Mn{M,S)dMd5 = pi, eq.(l) 
gives the physical number density of region (M, 6) at 
ti- The comoving number density is n{M,S)com = 
n{M,5){ai/ao)^, where a{t) is the cosmic factor. 

It has been shown (Carr 1975) that when a region 
(M, 6) enters the horizon at th, i.e. when M is about 
the same as horizon mass at th, this region will col- 
lapse to a black hole if the initial density contrast S 
satisfies 
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where parameter 7 ~ 1/3 and Mj is the horizon 
mass at ti, i.e. the horizon mass at the beginning 
of the radiation-dominated. The mass of a black 
hole formed from this region is equal to Mi,h = 
7V2mV^M2/3. Thus, the mass function of PBHs is 
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where 6h = 5{M/Mi)'^/^ is the density contrast when 
the region (M, 8) crosses the horizon at th, and vari- 
ance Uh = aR{M/Mifl^. Since M > Mj, the mini- 
mum mass of the PBH is 

For a power law initial density perturbation with 
spectral index n, i.e. ah cx: we have 

51n(j/?/51nM = (n-|-3)/6, and eq.(3) becomes (Kim 
& Lee 1996) 



nbh{Mbh 
n + 3 /2\i/2 



(4) 
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As usual, the variance ah will be normalized with the 
COBE observation (e.g.. Fang & Xu 1999, Bugaev & 
Konishchcv 2001). 

In eq.(4), only the factors p., = 3/327rGtf and 
Mi = (47r/3)(cti)'^pi are dependent on tj. Thus, 

nbh{Mbh) oc t^^^^. On the other hand, {ai/aof ~ 



(1 + Zeqf'{ti/teqf'^'^ , where teq and Zeq are, respec- 
tively, the time and the redshift at the end of the 
radiation-dominated era. Thus, the comoving mass 
function {ai/ao)^nbh{Mbh) is independent of ti. 

2.2. Ionizing photon background from PBH's 
evaporation 

For a Schwarzschild black hole m, the Hawking 
emission rate for particles with spin s, energy from E 
to E + dE and per degree of particle freedom is 



dN 
'dt 



V,dE 



2lTn eE/kBT{m) _ (-l)2s' 



(5) 



where T{m) = hc^/SnGkB m is the temperature of 
black hole m, and the dimensionless parameter F^ is 
the absorption probability of spin s particle by a black 
hole m (MacGibbon k Webber 1990; Page 1975). 
For low energy photon emission, i.e., hv <C kbT(m), 
r = 128G^m^z/^/3ci2, and eq.(5) gives 



dvdt 



f{iy, m) 



3cl2[g/»i//fcsT(m) _ ij- 



(6) 



The evaporation leads to the decreasing of black 
hole mass with the rate 



dm 
It 



a{m) 



(7) 



where a{m) takes into account all degrees of freedom 
of evaporated particles (MacGibbon 1991). Since 
a(m) is not a strong function of m, the time depen- 
dence of m for a black hole with initial mass Mbh at 
th is approximated by 



miMbh,t)c^[Mgh-Mm)t]'/^ 



(8) 



if t > th. 

Thus, the background photon energy flux J{u, z^^bs) 
per unit frequency interval and per solid angle at red- 
shift Zobs is as follows: 



J {v, Zobs) _ _c_ 
hv Ait 



rt{zobs) { 

J tea V 



dt{^)(^) (9) 
a J \aobsJ 



/ dMbhUbhiMbh) exp[-re//(z^, Zobs, z)\ 
f[iy{l + z)/{l + Zobs),m{Mbh,t)], 

Mmin is given by max[7-^/^Mj, m(t)], where m{t) is 
the solution of m{t) = {3a[m{t)]t}^^'^ . The integral 
upper limit M^ax in eq.(9) in not important, as the 
mass function nbh{Mbh) is rapidly decreasing with the 
increasing of Mbh- Tej f{i^^ Zobs^ z) is the opacity for 
photons emitted at z and observed at Zobs at the fre- 
quency V. 
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3. RESULTS 

Using eq.(9), we calculate the ionizing photon back- 
ground at redshift Zobs ~ 6 - 8. The opacity Teff for 
photons observed at the hydrogen Lyman edge {hi^Q 
=13.6 ev) is dominated by HI photoelectric absorp- 
tion. We will use the Teff estimated by the so-called 
high absorption model (Meiksin Sz Madau, 1993). It 
is 

T{l^O,Zobs,z) (10) 

0.0097(1 + Zobs)'-''[{l + ^)'-"' - (1 + Zobs)'-'"] 

-0.0068(1 + Zo6s)'-''[(l + zf^ - (1 + Zobsf^] 

8mxlO-'[il + zf-^-il + Zobsf% 

The optical depth given by eq.(lO) is very short. For 
Zobs = 6, the mean transmission e~'^''ff is less than 
10~^ at z ~ 6.5. That is, the effective flux J{uq) 
observed at 2; = 6 is only contributed by PBHs local- 
ized within the spatial range from 2; = 6 to 6.5. For 
2; = 8, this range is about z = 8 to 8.25. Actually, 
eq.(lO) is found from a fitting with observations at 
z ~ 3 - 4. However, there is no direct observation 
of the hydrogen Lyman edge absorption at redshifts 
z > Q. Nevertheless, eq.(lO) is consistent with the 
Lya transmission = 0.004 zt 0.003 in the redshift 
range [5.95 - 6.15] (Becker et al. 2001). Hence, it 
would be reasonable to estimate the optical depth at 
z = 6 - 6.5 by eq.(lO). 

The latest observation of the complete Gunn- 
Peterson troughs in the absorption spectrum of a 
z=6.28 quasar indicates that the photoionization 
rate at z ^6 is lower than 8.0xl0~-'^^s~^ for Lya, 
or 2.0xl0-^^s-^ for Ly7 (Fan et al. 2001). As- 
suming the ionizing photons are from massive stars 
J(z^) oc i>~^ (Barkana & Loeb 2001), one can then 
find that the upper limits to the ionizing photon back- 
ground J-21 = J(i^)/10^^^ {Ji^) ill vldH of erg cm^^ 
sec~^Hz~^sr~^) at photon energy 13.6ev are, respec- 
tively, 5.3 X 10^^ and 1.3 x 10~^ for Lya, and Ly7, 
and at 24ev, are 3.1 x 10~^, 7.8 x 10~^ for Lya, and 
Ly7, respectively. 

Figure 1 presents the PBH contributed ionizing 
photon background J-21 at z = 6 as a function of 
the index n of the primordial fluctuation power spec- 
trum. The upper limits to the ionizing photon back- 
ground yield n < 1.273 for £'=13.6ev, and n < 1.270 
for i5;=24 ev. 

It should be pointed out that these constraint basi- 
cally are independent of U. Eq.(9) shows that J(fo) 
depends on tj only via M„iin, the truncation of PBH 
mass function at small mass side. On the other hand, 
the optical depth at 2; = 6 is very small. Only 
PBHs, which can survive to z ~ 6, has contribu- 
tion to J{i'o). Thus, J(fo) is independent of Mmim 
if Mjnin < 10^^ g, as those PBHs have already evap- 



orated at z ~ 6 [eqs.(7) and (8)]. In inflationary 

universe, U = 0.301^7^^^ Mp;/r4 (Kolb & Turner 
1990). Thus, the M^j„-independence leads to also 
Tr/i-indep endence . 

The Mmin or T^/i-indep endence of J(i'o) can be also 
seen from Fig. 2, which plots dJ{i'o, Zobs)/dm vs. m. 
It shows that J(i^o) at z = 6 is mainly contributed 
by the evaporation of PBHs with mass m ~ 2 x 10^*^ 
g. Therefore, J(fo) is independent of Mmin if it is 
less than 2 x 10-*^^ g. The peak 2 x 10^^ g. is also 
consistent with the assumption of hu -C ki,T{m) used 
in eq.(6). 

We calculate the flux J(z^, z) at z = 6 in the fre- 
quency range from 1/ = 3x10^^ to 6 X 10^^ which 
corresponds to photon energy range from 13.6 to 24.6 
ev (Hel ionization). In this frequency interval one 
can still approximately use the optical depth Tg// of 
eq.(lO). The result is shown in Fig. 3. The flux J(j^) 
can be well approximated by a power law spectrum 
J(jv) oc u^. This spectrum is very different from the 
models of massive stars (J(i^) oc v^^) and quasars 
(J(z^) oc see Madau, Haardt, & Rees 1999), for 

both of which, the flux decreases with the increas- 
ing of v. This is because J(j^) are mainly from the 
low energy tail of the photon emission by lO^^g black 
holes, whose temperature T(m) ^ huo/kb. In the low 
energy tail, the intensity is increasing with u. 

In all the above calculations, the photons produced 
from the fragmentation of the PBH-emitted quarks 
and gluons are not considered (MacGibbon 1991). 
Since J(fo) is given by local PBHs, we need only 
to consider the fragmentation of quarks and gluons 
into low energy (~ 10 - 20 ev) photons. This prob- 
lem is still quite unclear and uncertain. Since the 
fragmentation effect will lead to an increase of J(i'o). 
Therefore, all above-mentioned upper limits on n will 
be held. 

4. DISCUSSIONS AND CONCLUSIONS 

In the standard inflationary model with a power- 
law power spectrum of the primeval density pertur- 
bations, we calculated the ionizing background flux 
of the PBHs. With the latest observed ionizing pho- 
ton background flux at the redshift of reionization, 
we find that the spectral index of the primeval den- 
sity perturbations should be n < 1.27, assuming that 
the background photon flux are emitted by massive 
stars. The difference is negligible if the background 
flux coming from quasars. If the parameter 7 is 
taken to be 0.6 (Niemeyer & Jedamzik 1998), we have 
n <1.29. We also calculate J(z^) at 2; = 8. The re- 
sult is almost the same as Fig. 3. However, there 
is no available data on the J{uq) at this redshift. If 
the reionization is really onset at 2: ~ 6 - 8, J(fo) at 
2 = 8 would be significantly less than that at 2; = 6, 
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one may have more tight constraint. All these con- 
straints are independent of the reheating temperature 
Trh, if Tj-h > W^^ Gev. The abundance of black hole 
given by this constraint is < 10"'*, which is much 
stronger than the critical density constraint ^Ibh < 1- 
The current observations of the angular power spec- 
trum of the temperature fluctuations of cosmic mi- 
crowave background radiation (CMB) are found to 
be able to fit with a n = 1 power spectrum on an- 
gular scales from tens of degrees to sub-degrees. The 
PBH's constraint n < 1.27 shows that primeval den- 
sity perturbations might still be scale invariant on 
scales much less than the direct observation of the 



CMB. Moreover, in the CMB data fitting, the pa- 
rameters of the spectral index n and the optical depth 
of the CMB photons to the last scattering surface Tq 
are degenerated (Stompor et al. 2001). That is, even 
when the fitted result is n fitting — 1; the original in- 
dex may be n > 1 if Tc > 0. On the other hand, the 
constraint on n given by the PBH is independent of 
the optical depth Tc. 

Thanks to Hongguang Bi for stimulating discus- 
sions. PH thanks M. Gibilisco for sending him her 
papers and helpful discussions. PH is supported by a 
Fellowship of the World Laboratory. 
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Fig. 1. — The PBH's ionizing background flux J_2i as a function of spectral index n of the primeval density perturbations at 
z—6: 1. photon energy E=13.6 ev (solid); 2. 24 ev (dashed). The observed upper limits to J-21 at E=13.6 and 24 ev are shown, 
respectively, by the shaded areas at top-right and bottom left. The upper and lower boundaries of the two areas refer to the limits 
given by Lya and Ly7 absorption troughs. These limits are estimated with the model of massive stars. 

Fig. 2. — The differential ionizing background flux dJ{iyo,'m) /dm vs. m, with the power spectrum index n=1.27. The solid, 
dotted, and dashed lines refer to cases of z=6, 8, and 10, respectively. 

Fig. 3. — The spectrum of ionizing photon produced by PBHs. The shaded area is upper limits to the ionizing photon spectrum 
given by Lya, and Ly7 in the massive star model. 
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